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ABSTRACT 


The  study  presents  a  modified  mathematical  expression 
for  the  wave-spectrum  function  in  the  Fourier  representa¬ 
tion  of  the  wave  pattern  of  a  ship  advancing  at  constant 
speed  in  calm  water.  This  new  expression  is  obtained  from 
the  well-known  usual  expression  via  several  applications 
of  Stokes'  theorem  for  combining  the  integrals  along  the 
top  waterline  and  over  the  hull  surface  of  the  ship.  The 
modified  expression  for  the  wave-spectrum  function  is 
considerably  better  suited  than  the  usual  expression  for 
accurate  numerical  evaluation,  notably  for  evaluating 
the  short  divergent  waves  of  interest  for  remote  sensing 
of  ship  wakes,  because  the  significant  numerical  cancella¬ 
tions  occurring  between  the  waterline  and  hull  integrals 
in  the  usual  expression  are  automatically  and  exactly 
accounted  for  in  the  modified  mathematical  expression,  as 
is  demonstrated  mathematically  and  confirmed  numerically. 
Whereas  the  values  of  both  the  velocity  potential  and  its 
gradient  at  the  hull  are  required  in  the  usual  expression 
for  the  wave-spectrum  function,  the  new  expression  only 
involves  the  tangential  velocity  at  the  hull,  not  the 
potential.  This  new  expression  thus  defines  the  wave- 
spectrum  function  in  terms  of  the  speed  and  the  size  of 
the  ship,  the  hull  form,  and  the  tangential  velocity  at 
the  mean  hull  surface. 


ADMINISTRATIVE  INFORMATION 

This  study  was  performed  at  the  David  Taylor  Research  Center  and  was  funded 
under  the  Applied  Hydrodynamics  Research  Program,  Program  Element  61153N,  Task  Area 
BR02301 51  and  Work  Unit  1542-109. 
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INTRODUCTION 


Near-field  potential-flow  calculations  about  ships  advancing  at  constant  speeds 
in  calm  water  are  routinely  required  for  evaluating  their  hydrodynamic  characteris¬ 
tics,  in  calm  water  and  in  waves,  and  for  determining  the  required  propulsion  and 
control  devices.  Calculations  of  far-field  ship  wave  patterns  are  also  important  in 
connection  with  wave-resistance  predictions  and  remote  sensing  of  ship  wakes.  In 
particular,  the  latter  practical  application  requires  the  ability  of  determining  the 
short  divergent  waves  in  the  wave  spectrum  having  wavelengths  between  5  cm  and  40  cm 
associated  with  Bragg  scattering  of  the  electromagnetic  waves  in  typical  SAR  systems 
used  in  remote  sensing  of  ship  wakes.  No  meaningful  predictions  of  such  short  waves 
can  be  obtained  on  the  basis  of  currently  available  numerical  methods.  More 
generally,  numerical  predictions  of  the  steady  wave  pattern  at  large  and  moderate 
distances  behind  a  ship  are  notoriously  difficult  and  unreliable,  as  was  recently 
made  clear  at  the  Workshop  on  Kelvin  Wake  Computations  [l].  Ship  wave-resistance 
calculations  are  also  known  to  be  difficult  and  unreliable. 

Alternative  numerical  methods  have  been  developed  for  evaluating  near-field  flow 
about  a  ship,  that  is,  flow  at  the  hull  surface  and  in  its  vicinity.  These  include 
finite-difference  methods,  e.g.  Coleman  [2]  and  Miyata  and  Nishimura  [3],  and  the 
more  widely  used  boundary  integral  equation  methods,  also  known  as  panel  methods. 

The  latter  methods  can  be  divided  into  two  main  groups,  according  to  the  Green 
function  that  is  used.  These  two  groups  of  methods  are  the  Rankine-source  method 
and  the  Neumann-Kelvin  method,  which  are  based  on  the  simple  Rankine  (free-space) 
fundamental  solution  and  the  more  complex  Green  function  satisfying  the  linearized 
free-surface  boundary  condition,  respectively. 

The  Rankine-source  method  was  initiated  by  Gadd  [4],  Dawson  [5]  and  Daube  [6], 
and  has  since  been  adopted  by  many  authors.  The  Neumann-Kelvin  approach  has  a  long 
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history.  A  survey  of  recent  numerical  predictions  obtained  by  a  number  of  authors 
on  the  basis  of  the  Neuraann-Kelvin  method  may  be  found  in  Baar  [7],  This  study  and 
that  by  Andrew,  Baar  and  Price  [8]  also  contain  extensive  comparisons  of  the  authors' 
own  Neumann-Kelvin  numerical  predictions  with  experimental  data.  An  approximate 
solution,  defined  explicitly  in  terras  of  the  value  of  the  Froude  number  and  the  hull 
shape,  to  the  Neumann-Kelvin  problem  was  proposed  in  Noblesse  [9],  This  slender- 
ship  approximation  was  recently  used  by  Scragg  et  al.  [10]  and  compared  to  both 
Neumann-Kelvin  predictions  and  experimental  data  in  [7]  and  [8]  and  to  experimental 
data  in  [ 1 ]  and  [ 11 ] . 

The  aforementioned  alternative  numerical  methods  for  predicting  flow  in  the 
vicinity  of  a  ship  are  not  all  directly  suitable  for  predicting  the  wave  pattern  of 
a  ship  at  large,  or  even  moderate,  distances.  More  precisely,  the  finite-difference 
method  and  the  Rankine-source  panel  method  require  truncating  the  flow  domain  at 
some  relatively-small  distance  away  from  the  ship  and  therefore  can  only  be  used  for 
near-f.'eld  fl<">w  calculations.  (However,  these  near-field  flow  predictions  can  be 
used  as  input  to  the  far-field  Neumann-Kelvin  flow  representation  considered  in  this 
study.)  On  the  other  hand,  the  Neumann-Kelvin  theoretical  framework  is  equally 
suitable  for  near-field  and  far-field  flow  predictions.  Indeed,  the  far-field 
Neumann-Kelvin  flow  representation  is  a  simplified  particular  case  of  the 
corresponding  near-field  representation. 

The  problem  considered  in  this  study  is  that  of  evaluating  the  steady  wave 
spectrum  and  the  wave  pattern  of  a  ship  at  moderate  and  large  distances  behind  it 
in  terras  of  the  near-field  flow  on  the  hull  surface.  The  near-field  flow  thus  is 
assumed  known  for  the  purpose  of  the  present  study,  which  is  concerned  with  the 
prediction  of  the  steady  wave  spectrum  and  the  wave  potential  behind  a  ship  stern 
within  the  Neuraann-Kelvin  theoretical  framework  as  was  just  noted. 
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This  theory  expresses  the  wave  potential  in  terms  of  a  Fourier  representation, 
as  is  well  known  and  is  specifically  indicated  by  Eq.  (20)  in  this  study.  The  wave- 
spectrum  (or  wave-amplitude)  function  in  this  Fourier  representation  is  defined  by 
the  sum  of  an  integral  along  the  mean  waterline  and  an  integral  over  the  mean 
wetted-hull  surface.  This  expression  for  the  wave -spectrum  function,  given  by  Eqs. 
(21 )— (23),  is  quite  ill  suited  for  accurate  numerical  evaluation  because  the 
waterline  integral  and  the  hull  integral  in  Eqs.  (22)  and  (23)  largely  cancel  out, 
as  is  shown  further  on  in  this  study  and  is  illustrated  in  Fig.  2.  Errors  in  the 
numerical  evaluation  of  the  waterline  and  hull  integrals  cause  imperfect  numerical 
cancellations  between  these  integrals  and  corresponding  large  errors  in  their  sum. 
This  fundamental  difficulty  was  recognized  in  [9]  and  in  Barnell  and  Noblesse  [ 1 2  J , 
where  attempts  to  remedy  it  were  presented.  However,  these  ad  hoc  approximate 
numerical  remedies,  based  upon  combining  the  waterline  integral  with  the 
contribution  to  the  hull  integral  stemming  from  the  upper  part  of  the  hull  surface 
are  not  satisfactory,  as  is  attested  by  the  fact  that  a  very  large  number  of  panels 
is  required  for  obtaining  reasonably  accurate  numerical  results  [ 1 2 ] • 

A  conceptually  simpler  and  numerically  more  effective  remedy  is  presented  in 
this  study,  in  which  a  modified  mathematical  expression  for  the  wave-spectrum 
function  is  obtained  via  several  applications  of  Stokes'  theorem  for  combining  the 
waterline  integral  and  the  hull  integral.  This  new  expression  for  the  wave-spectrum 
function  is  considerably  better  suited  than  the  usual  expression  for  accurate 
numerical  evaluation,  notably  for  evaluating  the  short  divergent  waves  of  interest 
for  remote  sensing  of  ship  wakes,  because  the  significant  numerical  cancellations 
occurring  between  the  waterline  and  hull  integrals  in  the  usual  expression  (see 
Fig.  2)  are  automatically  and  exactly  accounted  for,  via  a  mathematical 
transformation,  in  the  modified  expression  obtained  in  this  study.  The  fundamental 
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advantage  of  the  new  expression  over  the  usual  one  is  apparent  from  Figs.  8a-e. 

Another  interesting  feature  of  the  modified  expression  for  the  wave-spectrum 
function  is  that  it  only  requires  the  tangential  velocity  at  the  hull,  not  the 
potential,  whereas  the  usual  expression  requires  the  values  of  both  the  velocity 
potential  and  its  gradient  at  the  hull.  The  modified  expression  thus  defines  the 
wave-spectrum  function  in  terras  of  the  speed  and  the  size  of  the  ship,  the  hull 
form,  and  the  tangential  velocity  at  the  mean  hull  surface.  This  expression  is 
suitable  for  use  in  conjunction  with  a  boundary-integral-equation  method  based  on  a 
source  distribution  or  any  other  numerical  method  in  which  the  velocity  vector  (but 
not  the  potential)  is  determined  on  the  mean  hull  surface.  It  provides  a  coupling 
between  a  far-field  Neumann-Kelvin  flow  representation  and  a  near-field  flow 
calculation  method  based  on  the  use  of  Rankine  sources  or  finite  differences,  in 
particular. 
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Fig.  1  Definition  Sketch 
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NEUMANN-KELVIN  REPRESENTATION  FOR  THE  STEADY  WAVE  POTENTIAL  OF  A  SHIP 


As  was  already  noted,  this  study  considers  steady  potential  flow  about  a  ship 
advancing  at  constant  speed  In  calm  water  of  Infinite  depth  and  lateral  extent. 
Nondlmenslonal  coordinates  and  flow  variables  are  defined  In  terms  of  the  length  L 
and  the  speed  U  of  the  ship  and  the  water  density  p.  The  undisturbed  sea  surface  is 
chosen  as  the  plane  z  =*  0,  with  the  z-axis  pointing  upwards,  and  the  x-axis  is  taken 
in  the  ship  centerplane  (port-  and  starboard-symmetry  is  assumed)  and  pointing 
towards  the  bow,  as  is  depicted  in  Fig.  1.  The  Froude  number  and  its  inverse  are 
denoted  by  F  and  v,  respectively,  and  are  given  by 

F  -  U/(gL)l/2  =  1/v  ,  (1) 

where  g  is  the  acceleration  of  gravity. 

Within  the  so-called  Neumann-Kelvin  theoretical  framework,  the  velocity  poten- 
tial  <J>(£),  at  any  point  C  “  (£  ,  n ,  C  0)  strictly  outside  the  ship  hull  surface,  is 


defined  by  the  following  integral  representation  [9]: 


♦  (f)  -  ♦(€  )  +  X(f;« ♦)  , 

(2) 

♦  -  F 2 J  G  n  2t  dJL  +  f  G  n  da  , 

J  w  x  y  J  h  x  ’ 

(3) 

X  -  F2/  [G(t  ♦  +s  ♦  )-<}>3G/3x]t  dl 
w  x  c  x  s  y 

-  <j)3G/3n  da  +  F2 /  G  tt  (<j>  )  dxdy  , 

(4) 

where 

it  (<J> )  *  [  3 4> / 3  x—  ( V <|>  )2 /2]  3  (3 4>  /3  z+F2 3 2 4>  /3  x2  )/3z 

-  3(V<J,)2/3x  +  V<J>-V (V<j>  )2 /2  +  0(fV  )  . 

(5) 

In  Eqs.  (3)  and  (4),  the  symbols  w,  h  and  f  represent  the  positive  halves  of  the 

mean  waterline,  of  the  mean  hull  surface  and  of  the  mean  free  surface,  respectively, 

as  is  depicted  in  Fig.  1  (where  h  *  s  +  b  with  s  ■  hull  side  and  b  =  hull  bottom). 

Furthermore,  di  is  the  differential  element  of  arc  length  of  w  and  da  the  differen- 

tial  element  of  area  of  h.  Also,  n*(n,n,n)is  the  unit  vector  normal  to  h 

x  y  z 
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and  pointing  outside  the  ship,  t  =  (t  ,  t  ,  t  =  0)  is  the  unit  vector  tangent  to  w 

x  y  z 
-4 

and  pointing  towards  the  bow  and  s  =  (s  ,  s  ,  s  )  is  a  unit  vector  tangent  to  h  and 

x  y  z 

pointing  downwards,  as  is  shown  in  Fig.  1.  In  Eq.  (4),  <p  and  <f>  represent  the 
components  of  the  velocity  vector  V 4>  in  the  directions  of  the  tangent  vectors  t  and 


s  to  h,  respectively. 

The  nonlinear  termn(<j>)  defined  by  Eq.  (5)  and  the  corresponding  free-surface 
integral  in  Eq.  (4)  are  associated  with  the  nonlinearities  in  the  free-surface 
boundary  condition.  This  nonlinear  term  is  usually  neglected  in  practice. 

—  —  -4  4 

The  term  G  =  G(£;x)  in  Eqs.  (3)  and  (4)  is  the  Green  function  for  port-  and 
starboard-symmetry  defined  as 

_  ___  4  4  4  -4 

G  =  G(C;x)  =  G(£ ;x,y , z )  +  G(£;x,-y,z)  ,  (6) 


-4  -4 

where  G(£;x)  is  the  Green  function  associated  with  the  linearized  free-surface 

2  2  2 

boundary  condition  3G/3 5  +  F  3  G/3£  =  0.  The  function  G(C;x)  represents  the 


linearized  flow  created  at  the  point  £  *  (K ,  n »  C  <  0)  by  a  unit  outflow  at  the 
point  x  =  (x,  y,  z  _<  0),  stemming  from  a  submerged  source  if  z  <  0  or  from  a  flux 
across  the  mean  free  surface  if  z  =  0  as  is  shown  in  Noblesse  [  1 3  ] .  In  the  fore¬ 
going  equations  and  hereafter,  £  represents  the  "calculation  point",  where  the 
potential  is  evaluated,  while  x  represents  the  "integration  point"  in  the  integrals 
on  w,  h  and  f. 

The  Green  function  may  be  expressed  as  the  sum  of  three  terms,  as  follows 
[13,9]: 

4ttG(£;x)  -  S(5;x)  +  2v2N(X)  +  4v2H(x-S)W(X)  ,  (7) 

-4 

where  v  is  the  inverse  of  the  Froude  number  given  by  Eq.  (1),  and  X  is  defined  as 

X  =  (X,  Y,  Z  <  0)  =  [v2(5-x),  v2(n-y),  v2(c+z)]  ,  (8) 

and  represents  the  vector  of  coordinates,  rendered  nondimensional  with  respect  to 
the  characteristic  wavelength  U  /g  instead  of  the  ship  length  L,  joining  the  free- 
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surface  mirror  image  (x,  y,  -z)  of  the  singular  point  (x,  y,  z)  to  the  calculation 
point  (£,  n,  ?)• 

The  terra  S(|;x)  in  Eq.  (7)  corresponds  to  a  superposition  of  a  Rankine  source  at 
the  singular  point  (x,  y,  z)  and  a  Rankine  sink  at  its  free-surface  mirror  image 
(x,  y,  -z),  as  follows: 

S(|;x)  =  -[  (£-x)2+(n-y)2+(?-z)2]  +[(£-x)2+(n-y)2+(?+z)2]  .  (9) 

The  second  terra  N(X)  in  Eq.  (7)  represents  a  nonoscillatory  near-field  (local)  flow 

disturbance.  Finally,  the  terra  W(X)  represents  the  system  of  Kelvin  waves  trailing 
behind  the  singular  point  (x,  y,  z),  as  is  indicated  explicitly  by  the  Heaviside 
unit-step  function  H(x-£). 

By  using  expression  (7)  for  the  Green  function  in  Eqs.  (6)  and  (3)-(4)  we  may 

express  the  potential  <}>  (£ )  defined  by  Eq.  (2)  as  the  sura  of  three  potentials,  as 
follows: 

<t>(f)  -  <t>s(!)  +  <|>N(I)  +  $„(!)  ,  (10) 

which  are  readily  defined  by  Eqs.  (2)-(4)  and  (6)  in  which  the  Green  function  G  is 

2  2 

simply  replaced  by  S/(4ir),  v  N/(2ir)  and  v  H(x-£)W/tt,  respectively. 

Numerical  evaluation  of  the  potential  4>  associated  with  the  singular  algebraic 
terra  S  given  by  Eq.(9)  is  a  fairly  simple  task  (especially  for  points  f  strictly 
outside  the  ship  since  the  terra  S  then  is  never  singular)  for  which  extensive  expe¬ 
rience  is  available  in  both  aerodynamics  and  hydrodynamics.  In  particular,  only  the 
hull-surface  integrals  in  Eqs.  (3)  and  (4)  need  be  considered  for  the  potential  <p 

O 

since  Eq.  (9)  shows  that  we  have  S  =  0  on  w  and  f,  where  z  =  0. 

Numerical  evaluation  of  the  nonoscillatory  near-field  potential  in  Eq.  (10) 
associated  with  the  nonoscillatory  near-field  terra  N  in  Eq.  (7)  is  also  a  relatively 
simple  task  because  the  terras  N(X)  and  VN(X)  are  sufficiently  well-behaved  functions 
that  can  be  evaluated  numerically  with  satisfactory  accuracy  and  efficiency.  In 
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*►  •¥ 

particular,  the  function  N(X)  remains  finite  at  the  origin  X  =  0  whereas  the  func- 
*► 

tlon  VN(X)  has  a  relatively  weak  singularity,  which  is  given  in  Noblesse  [14], 
Furthermore,  this  singularity  does  not  occur  for  points  £  strictly  outside  the  ship. 
The  integrals  on  w,  h,  and  f  in  Eqs.  (3)  and  (4),  where  G  is  replaced  by  N,  can  then 
be  evaluated  numerically  by  using  ordinary  Integration  rules  since  the  integrands  in 
these  integrals  are  continuous  everywhere  and  nonoscillatory  (except  for  the  poten¬ 
tial  <p  and  the  free-surface  nonlinear  term  7r ,  which  are  wavy  functions  with  a  char¬ 
acteristic  wavelength  explicitly  defined  in  terms  of  the  value  of  the  Froude 
number).  Accurate  and  efficient  methods  for  numerically  evaluating  the  terms  N  and 
VN  have  recently  been  developed  by  Newman  [ 1 5 ]  and  by  Telste  and  Noblesse  [ 1 6 ] . 
Numerical  evaluation  of  the  wave  potential  $  in  Eq.  (10)  associated  with  the 

wave  terra  W  in  Eq.  (7)  is  considered  in  this  study.  The  wave  potential  *  ,  dominates 

w 

the  nonoscillatory  near-field  potential  4  +  <j>  at  some  distance  behind  the  ship,  as 

o  In 

is  well  known,  and  <j>w  is  the  most  important  component  for  practical  applications  to 
wave-resistance  predictions  and  remote-sensing  of  ship  wakes.  The  behavior  of  the 
oscillatory  term  W,  which  represents  the  system  of  Kelvin  waves  trailing  behind  the 

singular  point  x  as  was  already  noted,  is  considerably  more  complex  than  that  of  the 
nonoscillatory  near-field  term  N.  In  particular,  the  functions  W(X)  and  VW(X)  have 
strong  and  complex  singularities  at  the  origin  X  =  0  and  are  ill  behaved  in  the 
vicinity  of  the  line  Y  »  0  =  Z  and  X  _<  0,  as  is  shown  in  Ursell  [ 1 7 ]  and  Euvrard 
[ 1 8 ] .  Accurate  numerical  evaluation  of  the  functions  W  and  VW  is  difficult  in  the 
vicinity  of  the  origin  and  more  generally  at  the  plane  X  =  0.  In  spite  of  these 
difficulties,  a  method  for  numerically  evaluating  the  wave  terms  W  and  VW  has 
recently  been  developed  by  Baar  and  Price  [  1 9  ]  ,  Ursell  [20]  and  Newman  [ 2 1  ] .  In 
principle,  this  method  could  be  used  for  evaluating  the  wave  potential  $w  in  a 
manner  analogous  to  that  briefly  explained  in  the  foregoing  for  evaluating  the 
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nonoscillatory  near-field  potential  $  .  However,  it  is  not  clear  a  priori  that  the 

previously-noted  complex  behavior  of  the  functions  W  and  VW  along  the  line  Y  *  0  =  Z 

and  X  _<  0  would  not  cause  serious  numerical  difficulties  for  evaluating  the  wave 

potential  at  the  mean  free-surface  plane  z  ■  0. 

An  alternative,  indirect  or  Fourier-type,  method  for  evaluating  a  is  used  in 

w 

[12]  and  this  study.  The  method  is  based  upon  the  following  Fourier  integral 


representation  [ 1 3 ]  of  the  wave  terra  W(X): 

W(X)  =  /  Im  exp[z(l+t2)+i(X+Yt)(l+t2)l/2]  dt  ,  (11) 

—CO 

where  Im  denotes  the  imaginary  part.  By  using  Eq.  (8)  in  Eq.  (11)  we  may  obtain 
W(X)  -  /  Im  E*(t ;0  E(t;x)  dt  ,  (12) 

—OO 

where  the  functions  E  (t;5)  and  E(t;x)  are  defined  as 

E*(t;U  -  exp[p2{c  +  i(u£  +  vn  )}  ]  ,  (13a) 

E(t;x)  -  exp[p2{ z  -  i(ux  +  vy)}]  ,  (13b) 

In  Eqs.  (13a, b)  we  have 

P  =•  vp  with  p  =  (l+t2)1^2  ,  (14a, b) 

u  ■  1/p  and  v  =  t/p  ;  (15a, b) 

we  then  have 

1  u  0  and  0  _<  v  _<  1  for  0  _<  t  £  00  ,  (16a, b) 

with  u2+v2*l.  (17) 


Equations  (12)  and  (13a, b)  yield 

W(5;x)  -  W(5;x,y,z)  +  W(C;x,-y,z)  ■ 

00 

2  2  2  2 
2  Jq  Im  exp(v  ?p  )  cos(v  npt)  exp(iv  £p) 

exp(P2z)  [ E+(t;x,y )  +  E_(t;x,y)j  dt  ,  (18) 

where  the  functions  E+(t;x,y)  are  defined  as 

R+(t;x,y)  -  exp[-iP2 (ux  ±  vy)]  .  (19) 

The  wave  potential  is  defined  by  Eqs.  (2)-(4)  where  the  Green  function  G 


2  — 

is  replaced  by  the  wave  terra  v  H(x-£)W/tt,  as  was  already  noted.  In  this  study,  we 
limit  our  attention  to  the  wave  pattern  behind  the  ship  stern;  we  then  have  x  >  £ 
and  H(x-£)  =  1  for  points  (x,  y,  z)  on  the  hull  surface  h  +  w.  The  indirect, 
Fourier-type  method  used  in  this  study  for  evaluating  the  wave  potential  consists  In 
interchanging  the  order  of  integration  between  the  integration  point  x  in  the 
integrals  in  Eqs.  (3)  and  (4)  and  the  Fourier  variable  t  in  the  integral  (18).  The 
wave  potential  <j>  is  then  expressed  in  the  form 

4>W(C )  ■  (2/ir )  JQ  exp(v2£p2)  cos(v2npt)  Im  exp(iv2£p)  [K+(t)  +K_(t)]  dt  .  (20) 


In  this  Fourier  representation,  the  functions  K+(t)  are  given  by  the  sum  of  the 


terms  ^  and  X  defined  by  Eqs.  (3)  and  (4),  where  the  Green  function  G  is  replaced  by 
2  2 

the  terms  v  exp(P  z)  E+(t;x,y). 


The  functions  K+(t)  may  then  be  expressed  in  the  form 

K(t)  =  K  ±(t)  +  K  ±(t)  , 

-  0  9 

K  ±(t)  =  /  E.  n  2t  d£  +  v2  f,  exp(F2z)E  n  da  , 

o  w  ±  x  y  JhH  ±  x 

R1±(t)  =  /  E.  (t  d>  +s  <t>  +iv2 p<})  )t  U 
<j>  1  w  ±  xr  t  xT  s  y 

+  v2P2  exp(P2  z  )E+  4>n+  da  +  /  ^  E+  ir  (<J> )  dxdy  , 

where  Eq.  (19)  was  used,  the  term  n+  is  defined  as 
m  =*  -n  +  i(un  ±  vn  )  , 

±7.  V  V  ' 


(21) 

(22) 

(23) 

(24) 


z  x  y 

and  f^  represents  the  portion  of  the  mean  free-surface  plane  upstream  from  the  plane 


x  =  C. 

The  wave  potential  <j>  (f )  in  Eq.  (20)  thus  Is  defined  in  terms  of  a  familiar 

w 

Fourier  superposition  of  elementary  plane  waves  propagating  at  angles  0  from  the 
x-axis  given  by 

tan  0  =  ±v/u  =  ±t  .  (25) 

The  amplitudes  of  these  elementary  plane-wave  components  are  essentially  given  by 
the  functions  K+(t),  which  may  thus  be  referred  to  as  the  far-field  wave-amplitude 
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functions  or  as  the  free-wave  spectrum  functions.  These  functions  contain  essential 
information  directly  relevant  to  a  ship's  wave  resistance  and  signature.  In 
particular,  the  wave  resistance,  R  say,  experienced  by  the  ship  is  defined  in  terms 
of  the  wave-spectrum  functions  K+(t)  by  means  of  the  well-known  Havelock  formula 

wR/(pU2L2)  =  /“  [K+(t)  +K_(t)]2p  dt  .  (26) 

The  important  information  contained  in  the  functions  K+(t)  indeed  represents  a 
significant  advantage  of  the  indirect  Fourier-type  method  over  the  direct  integra¬ 
tion  method  mentioned  previously.  Furthermore,  the  integral  (20)  defining  the  wave 
potential  4>w(5)  is  more  amenable  to  numerical  evaluation  than  the  integral  (11) 
defining  the  wave  term  W(X)  because  the  wave-spectrum  function  K+(t )  +  K_(t) 
vanishes  as  t  -*•  «. 

Numerical  evaluation  of  the  wave  integral  (20)  and  of  the  functions  Kg~(t)  and 

+ 

K  (t)  defined  by  Eqs.  (22)  and  (23)  are  the  two  main  numerical  tasks  which  must  be 
9 

considered  in  the  indirect  Fourier-type  approach.  The  second  of  these  tasks  is 
examined  here.  Numerical  evaluation  of  the  waterline,  hull  and  free-surface 
integrals  in  Eqs.  (22)  and  (23)  is  a  seemingly  relatively-simple  task,  given  the 
value  of  the  potential  $  on  the  mean  hull  surface  h  and  waterline  w  and  the  value  of 
the  free-surface  nonlinear  term  tt  (<(>);  in  particular,  the  integrands  of  the  integrals 
in  Eqs.  (22)  and  (23)  are  continuous  functions. 


Nevertheless,  accurate  and  efficient  numerical  evaluation  of  the  functions 
+  + 

(t)  and  K,  (t)  requires  careful  analysis  because  the  trigonometric  functions 

U  (p 

2  2  2 

E+(x,y;t)  defined  by  Eq.  (19)  oscillate  very  rapidly  for  large  values  of  P  =  v  p  , 


as  is  the  case  for  typical  values  of  the  Froude  number  F  =  1/v  and  of  the  Fourier 
2  2  2 

variable  p  ®  1+t  =  sec  0,  and  because  the  potential  <j>  in  the  integrands  of  the 

2 

waterline  and  hull  integrals  in  Eq.  (23)  is  multiplied  by  the  large  numbers  v  p  and 
2  2  2  2 

v  P  *  (v  p)  ,  respectively.  The  terras  involving  the  potential  <j>  in  F,q.  (23)  there- 
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2 

fore  are  dominant  and  quite  large  for  typical  values  of  v  p.  More  precisely,  let 


Eq.  (23)  be  expresed  in  the  form 

%  -  kw  +  +  *V =  kw  +  io (V "  *  (27) 

where  the  free-surface  integral  and  the  superscript  ±  have  been  ignored  for 

simplicity  and  the  functions  K,,  ,  K and  K are  defined  as 

W  W  H 

Kw  ■  /„  E±  (tA+VS)Cy  “  •  <28»> 

V  "  /w  E±  +Cy  84  »  (28b) 

Kh"  *  fh  exp(P2z)E+  4>n±  da  ,  (28c) 

with  a  defined  by 

o  =*  v2p  *  (1+t2)1^2/  F2  ”  sec0  /  F2  .  (29) 


The  real  and  imaginary  parts  of  the  sum  of  the  port  and  starboard  contributions  to 
the  functions  ,  ioK^  and  o  K^-*  are  depicted  in  Fig.  2  for  0  _<  t  *  tan0  jC  10 

(corresponding  to  0  <  0  <  84°)  for  a  very  simple  case  corresponding  to  a  simple 
mathematicarry-def ined  hull  form  with  an  assumed  simple  mathematical  expression  for 
the  velocity  potential  <ji  at  the  hull.  More  precisely,  the  mathematical  hull  consid¬ 
ered  in  Fig.  2  has  constant  draft  and  rectangular  framelines,  with  draft/length  and 
beam/length  ratios  equal  to  0.07  and  0.16,  respectively.  The  hull  consists  of  a 
pointed  bow  region  0.2  _<  x  _<  0.5  with  parabolic  waterlines,  a  straight  middle-body 
region  -0.3  x  _<  0.2  and  a  rounded  stern  region  -0.5  <  x  j<  -0.3  with  elliptic  water- 
lines.  The  potential  in  Eqs.  (28a, b,c)  is  taken  as  <p  =  F  exp(v  z)  cos[v  (x-l/2)-3n/8], 
which  corresponds  to  an  elementary  plane  wave.  The  foregoing  simple  hull  form  and 
assumed  simple  expression  for  the  potential  at  the  hull  are  used  for  the  calcula¬ 
tions  presented  in  Fig.  2,  and  for  all  the  calculations  presented  further  on, 
because  this  simple  case  is  adequate  for  the  present  purpose  of  numerically  illus¬ 
trating  the  essential  properties  of  the  alternative  mathematical  expressions  for  the 
wave-spectrum  function  examined  in  this  study  and  it  permits  accurate  calculations 
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(the  required  integrations  can  be  partially  performed  analytically).  The  results 

presented  in  Fig.  2  correspond  to  a  value  of  the  Froude  number  equal  to  0.15. 

Figure  2  shows  that  the  function  K  is  considerably  smaller  than  its  components 

9 

Kw  ,  ioK^'  and  .  In  particular,  the  waterline  and  hull  integrals  laK^  and 

2 

a  K  '  do  not  appear  to  vanish  in  the  limit  9  -*■  90°  and  the  waterline  integral 
ri  W 

vanishes  appreciably  more  slowly  than  the  function  K  .  Significant  cancellations 

9 

therefore  occur  between  the  functions  iaKw'  and  and  between  the  sum  of  these 

two  functions  and  the  function  K^.  These  significant  cancellations  occur  for  all 

values  of  9  but  are  especially  notable  for  large  values  of  9,  which  correspond  to 

the  short  divergent  waves  in  the  spectrum.  The  errors  which  inevitably  occur  in  the 

numerical  evaluation  of  the  components  ,  lcrK^'  and  cause  imperfect  numerical 

cancellations  between  these  components  and  corresponding  large  errors  in  their  sum. 

Numerical  errors  in  the  sum  K  can  be  especially  difficult  to  control  because  the 

9 

errors  associated  with  the  numerical  evaluation  of  the  hull  integral  "  and  the 

H 

waterline  integral  +  laKy””  are  not  necessarily  comparable  (due  to  differences  In 

the  errors  associated  with  numerical  integration  over  hull  panels  and  waterline 

segments).  The  usual  expression  (23)  for  the  Neumann-Kelvin  correction  K  in 

9 

Eq.  (21)  thus  is  quite  ill  suited  for  accurate  numerical  evaluation.  A  modified 
mathematical  expression  in  which  the  cancellations  between  the  waterline  and  hull 
integrals  depicted  in  Fig.  2  are  automatically  and  exactly  accounted  for  is 
presented  further  on  in  this  study. 
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MODIFIED  EXPRESSIONS  FOR  THE  FAR-FIELD  WAVE-SPECTRUM  FUNCTION 


Modified  expressions  for  the  functions  K*(t)  and  K  *(t)  can  be  obtained  by 

U  <p 

using  Stokes'  theorem 

fc  V*t  dl  -  /  (VxV)  •  n  da  ,  (30) 

where  V  represents  a  vector  field,  and  t  and  n  are  unit  vectors  tangent  to  a  closed 


curve  C  and  normal  to  an  open  surface  S  bounded  by  C,  respectively.  In  particular, 


we  will  use  two  special  forms  of  Stokes'  theorem  corresponding  to  the  vector  fields 

V  ■  f  e  and  V  ■  f  e  ,  where  e  and  e  represent  unit  vectors  along  the  y-  and  x-axes 
y  z  y  z 

and  f  stands  for  a  scalar  function,  namely, 

/_  t  f  dl  ■  /_  (n  3f/3x  -  n  3f/3z)  da  ,  (31a) 

4  c  y  o  z  x 

f  t  f  dt  ■  /  (n  3f/3y  -  n  3f/3x)  da  .  (31b) 

u  z  b  x  y 

The  functions  K^^t),  which  correspond  to  the  so-called  zeroth-order  slender-ship 
approximation  to  the  Neumann-Kelvin  theory  [9],  are  considered  first,  and  a 
numerically-convenient  modified  form  of  Eq.  (22)  is  obtained. 
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The  zeroth-order  slender-shlp  approximation 


It  is  convenient  to  divide  the  mean  hull  surface  h  into  two  parts,  namely  the 
hull  side,  s  say,  and  the  hull  bottom,  b,  as  is  depicted  in  Fig.  1.  Equation  (22) 
then  becomes 

K  "(t)  =  /  E  n  2t  d£  +  v2/  exp(P2z)E.  n  da 
0  Jw±xy  Js  ±x 

+  v2/,  exp(P2z)E  n  da  .  (32) 

b  i  x 

The  hull  bottom  of  a  typical  ship  is  a  nearly  horizontal  surface,  so  that  we  have 
n^  -  0  on  b,  but  n^  is  usually  significant  on  the  hull  side  in  the  bow  and  stern 
regions.  However,  the  hull  side  of  a  typical  ship  is  a  nearly  vertical  surface, 
i.e.  we  have  n^  =  0  on  s.  It  is  therefore  convenient  to  express  the  integral  on  the 
hull  side  in  Eq.  (32)  as  an  integral  Involving  the  source  density  n^  by  using 
Stokes'  theorem  in  the  form  of  Eq.  (31a). 

More  precisely,  Eq.  (31a),  in  which  the  open  surface  S  and  the  function  f  are 

2  2  2 

taken  as  the  hull  side  s  and  the  functions  F  u  exp(P  z)E+,  yields 

/  exp(P  z)E  n  da  =  -iu  /  exp(P  z)E+  n  da 

-  F2u2  J  E+  t  d£  +  F 2u2  J  exp(P2z)E+  t  d£  ,  (33) 

where  Eqs.(l),  (14a)  and  (15a)  were  used,  and  w'  is  the  waterline-like  curve 
separating  the  hull  side  and  the  hull  bottom,  as  is  shown  in  Fig.  1.  The  unit  tan¬ 
gent  vector  t  =  (t^,  ty,  t^)  to  w'  is  pointing  towards  the  bow.  In  Eq.  (33),  the 
identity  t^  =  0  along  the  stem  and  stern  lines,  which  lie  in  the  ship  centerplane 
y  *  0,  was  used. 

By  substituting  Eq.  (33)  into  Eq.  (32)  we  may  obtain 

K  *(t)  -  /  E  (n  2-u2)t  di  +  u2  /  ,  exp(P2z)E  t  d£ 

0  ;w±x  y  ±y 

-  iv2u  /  exp(P2z)E  n  da  +  v2  /  exp(P2z)E+  n  da  .  (34) 

Comparison  of  Eqs.  (32)  and  (34)  shows  that,  on  the  hull  side,  the  source  density  n^ 
In  Eq.  (32)  has  been  replaced  by  the  density  ~iunz  Ecl*  (34).  The  latter  density 
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is  null  for  a  wall-sided  ship  and,  more  generally,  vanishes  in  the  limit  t  ♦  «,  as 
may  be  seen  from  Eqs.  (15a)  and  (14b).  The  hull-side  integral  therefore  is  gener¬ 
ally  less  important  in  the  modified  expression  (34)  than  in  the  original  expression 
(32).  In  particular,  the  hull-side  and  hull-bottom  integrals  in  Eq.  (34)  are  null 
for  a  wall-sided  ship  with  a  flat  horizontal  bottom,  for  which  Eq.  (34)  expresses 
the  functions  K^^t)  as  the  sum  of  two  line  integrals. 

The  trigonometric  functions  E+  defined  by  Eq.  (19)  are  rapidly  oscillatory,  as 
was  already  noted.  The  dominant  contribution  to  the  waterline  integral  in  Eq.  (34) 
therefore  stems  from  the  point(s),  if  any,  where  the  phases  P  (ux  ±  vy)  of  the  func¬ 
tions  E+  are  stationary.  These  points  of  stationary  phase  are  defined  by  the 
conditions  udx  ±  vdy  *  0,  which  yield  the  relations 

ut  ±  vt  -  0  ,  c  -  v  ,  t  «  +u  ;  (35a, b,c) 

*  y  a  y 

the  latter  two  relations  can  be  obtained  from  Eq.  (35a)  by  using  Eq.  (17)  and  the 

identity  t  2+t  2  *  1. 
x  y 

2 

The  terra  u  in  the  integrand  of  the  modified  waterline  integral  in  Eq.  (34) 

stems  from  the  hull-side  integral  In  Eq.  (32).  We  have  n  *■  -t  along  the  top  water- 

x  y 

2  2 

line  for  a  wall-sided  ship.  It  may  then  be  seen  from  Eq.  (35c)  that  the  terra  -u 
in  the  integrand  of  the  waterline  integral  in  Eq.  (34)  vanishes  at  a  point  of 
stationary  phase  for  a  wall-sided  ship.  This  result  indicates  that  the  waterline 
Integral  and  the  hull-side  integral  In  the  original  expression  (32)  cancel  one 
another  in  a  first  approximation  (more  precisely,  within  the  stationary-phase 
approximation)  for  a  wall-sided  ship.  The  major  contributions  stemming  from  these 
two  integrals  thus  are  combined  into  the  modified  waterline  integral  in  the  alterna¬ 
tive  expression  (34);  and  the  modified  hull-side  integral  in  Eq.  (34)  is  less  impor¬ 
tant  than  the  original  hull-side  integral  in  Eq.  (32),  as  was  already  noted. 

For  the  simple  mathematical  hull  form  defined  previously  we  have  n^  ■  0  on  the 
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hull  bottom  b  and  =  0  on  the  hull  side  s  and  Eqs.  (32)  and  (34)  then  yield 

=  K  +  K  ,  (36a) 

0  w  s 

K  *  K  *  +  K  ..  ,  (36b) 

0  w  w 

where  and  represent  the  waterline  and  hull-side  integrals,  respectively,  in 

the  usual  expression  (32)  and  K^*  and  correspond  to  the  integrals  along  the  top 

and  bottom  waterlines  w  and  w' ,  respectively,  in  the  modified  expression  (34).  The 

real  and  imaginary  parts  of  the  sum  of  the  port  and  starboard  contributions  to  the 

functions  K  ,  K  ,  K  *  and  K  >  are  depicted  in  Figs.  3a-e  for  0  <  tan0  <  5  (i.e.  for 
w  s  w  w  —  — 

0  _<  0  79°)  and  for  the  simple  hull  form  considered  previously  and  five  values  of 


the  Froude  number,  namely  F  =  0.1,  0.15,  0.2,  0.25  and  0.3.  For  F  =  0.1  and  0.15, 


Figs.  3a, b  show  that  the  bottom-waterline  integral  K  ^  is  quite  small  for  all  values 

w 

of  0  and  the  function  is  well  approximated  by  the  modified  top-waterline  integral 

K  *,  that  is  we  have  K  ^  «  K  *  =  K.„.  These  figures  also  show  that  the  modified 
w  w  w  0 

waterline  integral  K  *  is  appreciably  smaller  than  the  waterline  and  hull  integrals 

w 

and  in  the  usual  expression  (36a),  in  accordance  with  the  previous  theoretical 
considerations.  For  larger  values  of  the  Froude  number,  Figs.  3c-e  show  that  the 


bottom-waterline  integral  is  significant  for  smaller  values  of  0  but  vanishes 

rapidly  (exponentially)  for  increasing  values  of  0,  so  that  we  have  K  =  K  *  for 

U  w 

sufficiently  large  values  of  0. 
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0  1  2  3  4  5/0  1  2  3  4  5 


tanQ-+-  tane-*- 


Fig.  3c  The  functions  K  ,  K  ,  K  *  and  K  -  for  F  ■  0.2  • 

w  s  w  w 
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0  1  2  3  4  5/0  1  2  3  A  5 

tane-+  tan0-* 


Fig.  3e  The  functions  K  ,  K  ,  K  *  and  K  ^  for  F  »  0.3  . 

w  s  w  w 


-25- 


The  Neumann-Kelvin  approximation:  first  transformation 


The  linearized  Neumann-Kelvin  approximation  to  the  far-field  wave-amplitude 

funtions  K+(t)  is  given  by  the  sum  of  the  zeroth-order  slender-ship  approximation 

+  + 

Kg  (t)  defined  by  Eqs.(22)  or  (34)  and  the  correction  term  K.^  (t)  defined  by  Eq.  (23), 

where  the  free-surface  integral  is  neglected  as  is  indicated  in  Eqs.  (27),  (28a-c) 

and  (29).  A  modified  form  of  Eq.  (23)  is  now  obtained.  Stokes'  theorem  (31a),  in 

which  the  surface  S  and  the  function  f  are  taken  as  the  hull  surface  h  and  the 
2 

functions  exp(P  z)E+<j>,  yields 

f  E.  d>t  d£  ■*  f.  fn  3  exp(P2z)E  (j>  /3x  -  n  3  exp(P2z)E  $  /3z]  da  ,  (37) 

where  we  used  the  identity  ty  =  0  along  the  stem  line,  the  keel  line  and  the  stern 
line  which,  together  with  the  top  waterline  w,  border  the  hull  surface  h.  By 
substituting  Eq.  (37)  into  Eq.  (23),  and  upon  using  Eqs.  (24),  (19)  and  (17),  we  may 
obtain 


K  ~  (t)  -  /  E  (t  <|>  +s  <(>  )t  dA 
<f>  Jw  ±  xYt  xTs  y 

+  iv2p  exp(P2 z )E+  [  nz3<J>/3x-nx3<j>/3z±v2pt(uny+vnx)<}>]  da  .  (38) 

Equation  (38)  can  be  modified  further  by  using  Stokes'  theorem  (31b),  where  S 

and  f  are  taken  as  h  and  exp(P2z)E+4> ,  respectively.  The  contribution  of  the  top 

waterline  w  is  null  because  we  have  t  =  0  along  w.  The  contribution  of  the  stem 

line,  the  keel  line  and  the  stern  line,  which  lie  in  the  ship  centerplane  y  =  0,  to 

the  sum  K  +(t)  +K  (t)  can  also  be  shown  to  be  null.  We  then  have 
4> 

P2  exp(P2z)E+  (un^+vnx)<f>  da  =  i  exp(P2z)E+  (nx3<}>/3y-n^3ij>/3x)  da  .  (39) 

By  substituting  Eq.  (39)  into  Eq.  (38),  and  upon  using  Eqs.  (14a)  and  (15b),  we  may 
obtain 


K.±(t)  =  /  (t  <p+s  <p  )t  dJL  +  iv2p  /,  exp(P2z)E  a  da  , 

4>  Jw±  xrtxrsy  'h  t  t 

where  the  amplitude-functions  a+  are  defined  as 

a+  =  nz3ij>/3x  -  nx3<))/3z  ±  iv(nx3(j>/3y  -  ny3<J>/3x)  . 


(40) 


(41) 
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Equation  (40)  may  be  expressed  in  the  form 

%=kw  +  1oKh’  (42) 

where  the  superscript  ±  was  ignored  for  simplicity,  K  represents  the  waterline 

w 

integral  defined  by  Eq.  (28a),  a  is  given  by  Eq.  (29)  and  K  is  the  hull  integral  in 
Eq,  (40),  that  is  we  have 

Kh  -  /h  exp(P2z)E±  a±  da  .  (43) 

Comparison  of  Eqs.  (27)  and  (42)  shows  that  we  have 

Kh  -  Kw'  -  iaKH'  .  (44) 

Thus,  the  waterline  integral  and  the  hull  integral  -ioK^'  have  been  combined  into 
the  modified  hull  integral  K  via  two  applications  of  Stokes'  theorem  in  the  form 

it 

given  by  Eqs.  (31a, b).  The  real  and  imaginary  parts  of  the  sum  of  the  port  and 
starboard  contributions  to  the  functions  ioK  " ,  and  ioK  are  depicted  in 

Figs.  4a-e  for  0  tan0  10  (i.e.  for  0  £  0  <  84°)  and  for  the  simple  hull  form  and 

the  assumed  simple  expression  for  the  potential  at  the  hull  considered  previously  in 
Fig.  2.  Figures  4a-e  correspond  to  the  following  five  values  of  the  Froude  number 
F:  0.1,  0.15,  0.2,  0.25  and  0.3.  These  figures  show  that  the  waterline  and  hull 
integrals  ioK^'  and  o  are  considerably  larger  than  the  modified  hull  Integral 
ioK  .  Although  the  latter  integral  is  identical  to  the  sum  of  the  integrals 
ioK  '  and  "  ,  it  clearly  is  preferable  to  evaluate  ioK.  directly  rather  than  the 

W  H  H. 

2 

sum  of  the  integrals  ioKy*  and  a  .  The  modified  expression  for  the  Neumann-Kelvin 

correction  term  K  given  by  Eqs.  (40)  and  (41)  or  (42)  therefore  represents  a  signi- 
ficant  improvement  in  comparison  with  the  usual  expression  given  by  Eqs.  (23)  or  (27). 

The  cancellations  between  the  waterline  integral  ioKw'  and  t^e  hull  integral 
a  K  '  depicted  in  Figs.  4a-e  can  easily  be  explained  mathematically  for  a  wall-sided 

rl 

2  j, 

ship  form.  For  large  values  of  P  ,  the  major  contribution  to  the  hull  Integral  K 

H 

stems  from  the  upper  part  of  the  hull  surface  in  the  vicinity  of  the  waterline. 
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where  we  have  n  =  -t  ,  a  =  t  and  n  =0  for  a  wall-sided  ship.  More  precisely, 
xy  yxz  r  J  ’ 

Eqs.  (28c)  and  (24)  yield 


~  -i  /  E  4> (ut  +vt  )  dl  (°  exp(P2z)  dz  as  P2  +  “  . 
H  Jw±  yx  J-oo 


We  then  have 


>  2 
ioK  /  ~  u  /  E,  d>  (ut  +vt  )  d£  as  P  *  »  , 

H  1  w  ±  T  y  x  ’ 

where  Eqs.  (29),  (14 a)  and  (15a)  were  used.  By  using  Eqs.  (28b)  and  (17)  we 


then  obtain 


K  /  -  ioK  "  ~  v  /  E  $(vt  ±ut  ) 
W  H  1  w  ±  y  x 


dj l  as  P  -*■  “  . 


For  large  values  of  P  ,  the  trigonometric  functions  E  defined  by  Eq.  (19)  are 

rapidly  oscillatory  and  the  dominant  contributions  to  the  waterline  integrals  K  '  , 

ioKjj'  and  K^-ioR  '  defined  by  Eqs.  (28b),  (45)  and  (46),  respectively,  therefore 

stem  from  the  point(s),  if  any,  of  stationary  phase  of  the  trigonometric  functions 

E+  defined  by  Eqs.  (35a, b,c).  At  such  a  point  of  stationary  phase  the  terms  t  , 

u(ut^+vtx)  and  v(vt^±utx)  in  the  integrands  of  the  waterline  integrals  (28b),  (45) 

and  (46)  take  the  values  +u,  +u  and  0.  respectively,  which  demonstrates  that  the 

waterline  integral  and  the  hull  integral  -ioK^'  cancel  out  in  a  first 

2 

approximation  in  the  limit  P  -►  00  for  a  wall-sided  ship. 

The  real  and  imaginary  parts  of  the  sum  of  the  port  and  starboard  contributions 

to  the  modified  waterline  and  hull  integrals  K  and  ioK  ,  respectively,  and  their 

W  H 

sum  K  are  depicted  in  Figs.  5a-e  for  the  cases  considered  previously  in  Figs.  4a-e 

and  Fig.  2.  It  may  be  seen  from  Figs.  5a-e  that  the  waterline  and  hull  integrals  Kt1 

W 

and  ioK  are  appreciably  larger  than  their  sum  K  ,  especially  for  large  values  of 

H  <j> 

tanQ.  Significant  cancellations  therefore  occur  between  the  waterline  and  hull 
integrals  in  Eq.  (40).  Further  modifications  of  the  expression  for  the  function 
K  defined  by  Eq.  (40)  are  then  desirable  for  numerical  calculations.  These 


modifications  are  now  presented. 
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The  Neumann-Kelvln  approximation:  second  transformation 

Let  a,  8,  y  represent  three  constants  and  e  ,  e  ,  e  three  unit  vectors  along  the 

x  y  z 

x— ,  y-,  z-axes,  respectively.  We  have 

n  •  Vx(ae  +f5e  +ye  )d>  *  a(n  3<b/3z-n  3 4> /3 y ) 
x  y  z  T  yT  z  T  ' 

+  B(n  3<J>/3x-n  3<j>/3z)  +  y(n  9<f>/3y-n  3<(>/3x)  .  (47) 

z  x  x  y 

Equation  (47)  shows  that  the  amplitude  functions  a +  defined  by  Eq.  (41)  may  be 
expresed  in  the  form 

a.  =  n  •  Vx(e  ±ive  )<b  .  (48) 

±  y  z  T 

2 

Let  the  functions  exp(P  z)E+,  where  E+  are  the  trigonometric  functions  defined 
by  Eq.  (19),  be  denoted  as  e+.  We  then  have 

e+  *  exp[p2{z-i(ux±vy )} ]  and  (49) 

(50) 


Ve.  =  -P  e. (iue  live  -e  )  . 

±  ±  x  y  z 

We  have 

e  Vx(ae  +8e  +ye  )<j>  *  Vxe+ (ae  +8e  -  Ve+x(ae  +8e  +ye  )$  - 

Equation  (50)  yields 

-Ve  x(oie  +8e  +ye  )d>  =  P2e.4m.  with 
+  X  y  2  x  x 

m  =  (8±ivY)e  -  (a+iuy)e  +  i(u8+va)e  . 

—  x  y  z 

Equation  (53)  shows  that  we  have 

ra+  =  0  if  a  *  -iuy  and  8  =  +ivy  , 
which  yields 

ae  +8e  +ye  =  -iy(ue  ±ve  +ie  )  . 
x  y  z  x  y  z 

->■  -►  + 

This  condition  merely  expresses  that  the  vector  aex  +  8e^  +  and  the  vector  Ve+ 


(51) 

(52) 

(53) 

(54) 

(55) 


defined  by  Eq.  (50)  are  colinear.  Equation  (51)  then  becomes 

e.Vx(ue  ±ve  +ie  )6  *  Vxe.  (ue  ±ve  +ie  )<b  . 

±  xyz  ±xyz 

Equations  (48)  and  (56)  yield 

e.a.  =  e,n  •  Vx(e  ±ive  )d>  , 

±  ±  ±  y  z 


(56) 


(57) 
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(58) 


e  n  •  Vx(ue  ±ve  +ie  )<j>  =  n  •  Vxe  (ue  ±ve  +ie  )d>  . 

±  x  y  z  ±  x  y  z  y 

Let  the  amplitude  functions  a+  be  expressed  in  the  form 

a  =  b  +  c 
±  ±  ±  ’ 

with  b+  defined  as 

b+  =  n  •  Vx  (ae^+8  )<J>  , 

where  a,  3,  y  are  constants.  Equations  (59),  (48),  and  (60)  yield 

c.  =  -n  •  Vx[ae  +(3“l)e  +(y+iv)e  ]<J>  . 

—  x  y  z 

Equation  (58)  now  shows  that  we  have 

e.c.  =  -T  n  •  Vxe. (ue  ±ve  +ie  )d> 

±±  ±  ±  x  y  zT 


(59) 


(60) 


(61) 


if  the  constants  a,  3>  Y  are  chosen  as 

a  -  ur±  ,  3  -  l±vr±  ,  Y  =  i(r±±v)  (62a, b.c) 

where  T+  is  some  arbitrary  constant.  Equations  (60),  (62a, b,c)  and  (47)  then  yield 


b+  =  ur+(n  3<j>/3z-n  3<j>/3y)  +  (l±vr+)(n  3(j>/3x-n  3(j>/3z) 
—  —  y  z  —  z  x 

+  i(T+±v) (n  3 <}> /3 y-n  3<J> /3 x)  . 

—  x  y 

Equations  (59)  and  (61)  show  that  we  have 

^  ^ 

e.a.  =  e.b.  -  r  n  •  Vxe  (ue  ±ve  +ie  )d)  . 

±±  ±  ±  ±  ±  x  y  z 


(63) 


(64) 


By  using  Stokes'  theorem  (30)  we  may  now  obtain 

fh  e±a±  da  =  /h  e±b±  da  -  T±  Jc  e± (utx±vty+itz )<(>  dZ  ,  (65) 

where  the  curve  c  consists  of  the  waterline  w  plus  the  bow-keel-stern  line,  which 
lies  in  the  centerplane  y  =  0.  The  hull  integral  in  Eq.  (40)  has  thus  been 
expressed  in  Eq.  (65)  as  the  sum  of  a  hull  integral  and  a  waterline  integral  involv¬ 
ing  the  constants  T+,  which  are  arbitrary  and  can  then  be  selected  at  will.  The 
identities  (64)  and  (65)  do  not  involve  the  term  P  $  appearing  in  Eqs.  (52)  and  (51), 
Along  the  bow-keel-stern  line  we  have  y  =  0  and  t  =  0  and  Eq.  (49)  yields 
e+(utx±vty+itz)  =  (ut^+it^)  exp[ P2 (z-iux)]  . 

Equations  (40)  and  (65)  then  show  that  the  contribution  of  the  bow-keel-stern  line 
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(66) 


to  the  sum  K  +  +  K.  can  be  rendered  null  if  the  constants  T,  are  taken  as 
<P  v  1 

r+  -  +r  . 

Equation  (65)  then  becomes 

[,  e,a.  da  “I,  e.b.  da  ±  T  /  E.  (ut  ±vt  )$  dt 
J  h  ±  z  *  n  i  ±  x  y 

since  we  have  t  ■  0  and  e+  *  E+  along  the  waterline,  where  z  *  0. 

Equations  (19),  (14a),  (15a)  and  (1)  yield 
E.  (ut  ±vt  )dt  -  iF2u2dEJ. 

We  then  have 

±J  E  (ut  ±vt  )<j»  dt  -  ±iF2u2[  (E.*)w  -(E.*)  _  ] 

J  w  ±  x  y  Y  L  ±Y  bow  ±r  sternJ 

+  iF2u2  J  E±(3<j»/3t)  d t  . 

•f 

The  bow  and  stern  contributions  to  the  function  K  +  K,  are  null  because  y  ■ 

9  9 

E+  *  E_  at  the  bow  and  the  stern.  Equations  (67)  and  (68)  then  yield 

/ h  e±a±  da  a  Jh  e±b±  da  +  1f2 «2r  Jw  E±(3*/3t)  dt  . 

By  substituting  Eq.  (69)  into  Eq.  (40)  we  may  obtain  the  following  alternative 

expression  for  the  Neuraann-Kelvin  correction  K  "  to  the  spectrum  function: 

9 

k  1  (t)  -  /  e  [  (t  <t>  +s  ♦  )t  ±ura^/at]  dt 

q>  '  ■'w  ±l  xTt  xTs  y  J 

+  iv2p  /h  exp(P2z)E±  b±  da  , 

where  Eqs.  (49),  (19),  (15a)  and  (1)  were  used. 

By  substituting  Eq.  (66)  into  Eq.  (63)  we  may  obtain 

pb±  “  ?A±  » 

where  the  amplitude  functions  A+  are  given  by 

A+  ■  T(n  3<{>/3z-n  3$/3y  ±  B(n  3$/3x-n  3$/3z)  +  iC(n  3<j>/3y-n  3q» /3x) 

—  y  z  z  x  x  y 

with  B  and  C  defined  as 

B  ■  p(vr-l)  and  C  ■  p(r-v)  . 

These  relations  yield 

T  ■  Cu  +  v  and  B  *  Cv  -  u  . 


(67) 

(68) 

0  and 

(69) 

(70) 

(71) 

(72) 

(73a, b) 


-41- 


By  using  Eqs.  (7L)  ,  (72)  and  (73a, b)  we  may  then  express  Eq.  (70)  in  the  form 

K  _(t)  =  J  E  [  (t  4>  +s  4)  )t  ±  u(Cu+v)36/3 tl  di 

+  iv2/^  exp(P2z)E+  A+  da  ,  (74) 

where  the  amplitude  functions  A+  are  given  by 

A.  =  (Cu+vMn  a<p/dz-n  ct<p/3y) 

±  y  z 

±  (Cv-u)(nz3ij)/3x-nx34i/9z)  +  iC(nx34> /dy-n^dtp/dx)  .  (75) 

The  constant  C  in  Eqs.  (74)  and  (75)  is  arbitrary  and  may  then  be  selected  at  will. 

More  precisely,  C  may  be  chosen  as  an  arbitrary  function  of  t.  These  equations  thus 

define  a  one-parameter  family  of  alternative  mathematically-equivalent  expresions 

+ 

for  the  functions  K  (t). 

The  amplitude  functions  in  the  integrands  of  the  hull  integrals  in  the 

alternative  expressions  (40)  and  (74)  are  given  by  pa+  and  +A+,  respectively.  It 

may  seem  from  Eqs.  (41),  (14b)  and  (15b)  that  we  have  a+  =  0(1)  and  pa+  =  0(t)  as 

t  ■*  ”,  whereas  Eqs.  (14b),  (15a,b)  and  (75)  show  that  A+  =  0(1)  as  t  «  for  any 

finite  value  of  the  limit  C(”)  of  the  arbitrary  function  C(t).  The  hull  integral  in 

Eq.  (74)  therefore  vanishes  more  rapidly  than  the  hull  integral  in  Eq.  (40)  as  t  +  “. 

This  result  implies  that  the  waterline  integral  in  Eq.  (74)  likewise  vanishes  faster 

than  the  waterline  integral  in  Eq.  (40).  Indeed,  the  term  t^t  <j>t  in  the  integrand 

2 

of  the  waterline  integral  in  Eq.  (40)  is  replaced  by  the  term  (t  t  ±uv±Cu  )<(>  in 

x  y  t 

2  2  2 
Eq.  (74).  We  have  t  t  ±uv±Cu  ~  t  t  ±uv  as  t  +  ®,  with  an  error  0(u  )  =  0(t~  ),  and 

x  y  x  y 

t^t^iuv  =  0  at  a  point  of  stationary  phase  of  the  trigonometric  function  E+  defined 
by  Eqs.  (35a, b,c).  The  cancellations  occurring  for  large  values  of  t  between  the 
waterline  and  hull  integrals  in  Eq.  (40),  as  is  depicted  in  Figs.  5a-e,  may  then  be 

expected  to  be  significantly  reduced  in  the  alternative  expression  (74). 

2 

For  large  values  of  t,  Eqs.  (14a, b)  show  that  we  have  P  »  1.  The  major  con¬ 
tribution  to  the  hull  integral  in  Eq.  (74)  therefore  stems  from  the  upper  part  of 
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the  hull  surface  In  the  vicinity  of  the  waterline.  For  a  wall-sided  hull  we  have 

n^  *  -t  ,  n^  =  t^  and  n^  =  0  on  h  in  the  vicinity  of  w,  and  Eq.  (75)  becomes 

A+  *  [vtx+ut^+C(utx±vt^)]3ij)/3z  -  iC3$/3t  ,  (76) 

where  the  identity  t  3<J) /b x  +  t  3<J>/3y  *  3<J)/3t  was  used.  The  hull  integral  in  Eq.  (74) 

x  y 

2 

can  be  approximated  by  a  waterline  integral  for  P  »  1.  Specifically,  we  have 

Ti v2/h  exp(P2z)E±  A±  da  ~  /  E±  dl  ,  (77) 

where  is  given  by  3  +iu  A+.  Equation  (76)  then  yields 

A^  *  +iu2[vtx+uty+C(utx±vty)]34>/3z  +  Cu23«j>/3t  .  (78) 

We  may  choose  sx  =  0.  We  then  have  $  ■  3 <J> /3 1  for  a  wall  sided  hull  and  the  amplitude 

4- 

function,  A^  say,  in  the  integrand  of  the  waterline  integral  in  Eq.  (74)  becomes 

A  ±  -  (t  t  ±uv±Cu2)3<b/3t  .  (79) 

w  x  y 

The  major  contributions  to  the  waterline  integrals  in  Eqs.  (74)  and  (77)  stem 
from  the  point(s)  of  stationary  phase  of  the  trigonometric  functions  E+  and  from  the 
end  points,  that  is  the  bow  and  the  stern.  At  a  point  of  stationary  phase,  Eqs. 


(35a, b,c)  hold  and  Eqs.  (78)  and  (79)  become 

A^~  *  +  iu23<|>/3z  +  Cu23ij)/3t  , 

A  -  =  ±  Cu23ij)/3t  , 
w 


(80a) 

(80b) 


where  Eq.  (17)  was  used.  These  equations  show  that  the  sum  of  the  waterline  and  hull 
amplitude  functions  at  a  point  of  stationary  phase  is  independent  of  the  constant  C, 
as  must  be  true  in  general  (i.e.  for  all  values  of  t,  including  the  limit  t  •*  <*> 
considered  here),  and  is  given  by 

A^~  +  A^-  *  +iu23<f>/3z  ,  (81) 

which  stems  from  the  hull  integral  in  Eq.  (74). 

+  2 

Equations  (78)  and  (80a)  show  that  the  amplitude  function  A^  is  0(u  )  at  the 
bow  and  the  stern  and  at  a  point  of  stationary  phase.  The  dominant  contribution  to 
the  hull  integral  in  Eqs.  (74)  and  (77)  therefore  stems  from  the  point(s)  of 
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3 

stationary  phase  and  is  0(Fu  ).  Equation  (80b)  shows  that  the  stationary-phase 
contribution  to  the  waterline  integral  in  Eq.  (74)  likewise  is  0(Fu  ).  The  bow  and 
stern  contributions  to  the  waterline  integral  in  Eq.  (74)  stem  from  the  term 

txty3<t>/3t  (82) 

in  Eq.  (79),  which  is  0(1).  The  bow  and  stern  contributions  to  the  waterline  inte- 
gral  in  Eq.  (74)  therefore  are  0(F  u  )  and  dominate  the  0(Fu  )  stationary-phase 
contributions  to  the  waterline  and  hull  integrals  in  the  limit  t  «>.  The  waterline 
integral  in  Eq.  (74)  thus  dominates  the  hull  integral  in  this  limit. 

However,  an  exception  to  this  general  rule  occurs  if  the  amplitude  function 
defined  by  Eq.  (82)  is  null  at  the  bow  and  the  stern.  Such  would  be  the  case  for 
cusped  ends,  for  which  we  have  t  =  0.  A  more  realistic  case  is  that  of  a  round- 

y 

ended  ship  form,  e.g.  an  oil  tanke>',  for  which  we  have  t  =  0  and  9<j>/9t  =  0  at  the 
bow  and  the  stern.  For  such  hull  forms,  the  dominant  contributions  to  the  waterline 
and  hull  integrals  in  Eq.  (74)  stem  from  the  point(s)  of  stationary  phase  and  are 

3 

0(Fu  ),  as  follows  from  Eqs.  (80a, b).  These  equations  show  that  significant  cancel¬ 
lations  might  then  occur  between  the  waterline  and  hull  integrals  in  Eq.  (74)  unless 
the  function  C(t)  vanishes  in  the  limit  t  °°.  Accordingly,  we  impose  that  the 
arbitrary  function  C(t)  in  Eqs.  (74)  and  (75)  satisfy  the  condition 

C  -*■  0  as  t  00  .  (83) 

An  obvious  choice  for  the  function  C(t)  satisfying  Eq.  (83)  is 

C  -  0  .  (84) 

+ 

The  corresponding  expressions  for  the  functions  K  are  readily  obtained  from  Eqs. 

9 


(74)  and  (75).  These  equations  may  be  expressed  in  the  form 


%  -  V  +kh  *  (85) 

where  the  superscript  ±  was  ignored  for  simplicity  and  the  functions  "  and  K 
correspond  to  the  waterline  and  hull  integrals  in  Eq.  (74),  respectively.  The  real 
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and  imaginary  parts  of  the  sum  of  the  port  and  starboard  contributions  to  the  func¬ 
tions  K  and  ioK„  defined  by  Eqs.  (40)-(43)  and  the  functions  K,/'  and  K,/'  defined 
w  K  W  H 

by  Eqs.  (74),  (75)  and  (84),  (85)  are  depicted  in  Figs.  6a-e  for  0  _<  tanQ  _<  10  (i.e. 

for  0  <  e  84°)  and  frtr  the  cases  considered  previously  in  Figs.  2,  4a-e  and  5a-e. 

It  may  be  seen  from  Figs.  6a-e  that  the  functions  and  vanish  more  rapidly 

than  the  functions  K  and  ioK  for  increasing  values  of  t  =  tan6,  in  accordance  with 

the  foregoing  theoretical  considerations.  The  cancellations  occurring  between  the 

waterline  and  hull  integrals  K  and  ioK  for  large  values  of  tan0  thus  are  signifi- 

w  ri 

cantly  reduced  in  the  alternative  expression  K  ,  which  is  therefore 

preferable  to  the  expression  K  +  ioK  for  large  values  of  tanQ.  However,  the 
functions  and  are  appreciably  larger  than  the  functions  Kw  and  ioK^  for 

small  values  tanQ,  especially  in  Figs.  6a,b,c  corresponding  to  small  values  of  the 
Froude  number,  and  significant  cancellations  thus  occur  between  the  waterline  and 
hull  integrals  and  f°r  s'03!!,  values  of  tanQ.  The  expression  Kw  + 

IoK  therefore  is  preferable  to  the  expression  K  "  +K  "  for  small  values  of  ta^Q, 
whereas  the  reverse  holds  for  large  values  of  tanQ. 

The  amplitude  functions  in  the  integrands  of  the  waterline  integrals  in  Eqs. 

(40)  and  (74)  are  nearly  identical  for  small  values  of  t  if  we  have 

v  +  Cu  «  1  as  t  0  .  (86) 

Equations  (41)  and  (75)  show  that  the  amplitude  functions  a+  and  +uA+  in  the 
integrands  of  the  hull  integrals  in  Eqs.  (40)  and  (74)  likewise  are  nearly  identical 
as  t  ♦  0  if  condition  (86)  and  the  condition  l-u(u-Cv)  «  1  hold.  By  using  Eq.  (17), 
we  may  express  the  latter  condition  in  the  form  v(v+Cu)  <<  1,  which  is  Identical  to 
condition  (86).  This  condition  therefore  ensures  that  the  waterline  and  hull 
integrals  in  Eq.  (74)  are  nearly  Identical  to  the  corresponding  integrals  in  Eq.  (40) 
in  the  limit  t  >  0. 
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The  large-  and  sraall-t  conditions  (83)  and  (86)  are  satisfied  If  the  arbitrary 


function  C(t)  is  selected  in  the  form 
C  =  -Auv  . 


(87) 


Condition  (86)  then  becomes  v(l-Au  )  «  1  as  t  +  0.  Equation  (17)  shows  that  we 
2  3 

have  v(l-\u  )  ~  v  as  t  +  0  if 


A  *  1  as  t  0  . 


(88) 


By  substituting  Eq.  (87)  into  Eqs.  (74)  and  (75)  we  may  obtain 

K  “(t)  =  f  E  a.  di  +  iv2  /,  exp(P2z)E  A.  da  , 

<J>  ■'w±±  ;  h  ±± 

where  the  amplitude  functions  a+  a^d  A+  are  given  by 

a+  =  (tx(t,t+sx!t>s)t;y  ±  uv(l-Au2  )3<j> /3  t , 

A+  =  v(l-Au2)(n  3<(>/3z-n  3<j> /3 y ) 

+  u(l+Av2)(nz3<J)/3x-nx3(j)/3z) 

-  iAuv(n  3 <}> /3 y— n  d<p/dx)  . 
x  y 

An  obvious  choice  for  the  function  A(t)  satisfying  condition  (88)  is 


(89) 


(90a) 


(90b) 


A  =  1  . 


(91) 


The  corresponding  expressions  for  the  functions  K  are  readily  obtained  from 

<P 

2  2 

Eqs.  (89)  and  (90a, b),  where  we  have  1-u  =  v  by  virtue  of  Eq.  (17).  These  equa¬ 

tions  may  be  expressed  in  the  form 

%-kw*+kh‘-  (92) 


if  if 

where  the  superscript  ±  was  ignored  for  simplicity  and  the  functions  Kw  and 
correspond  to  the  waterline  and  hull  integrals  in  Eqs.  (89),  respectively.  The  real 
and  imaginary  parts  of  the  sum  of  the  port  and  starboard  contributions  to  the  func¬ 
tions  Ky  and  ioK^  defined  by  Eqs.  (40)-(43)  and  the  functions  K^*  and  defined  by 
Eqs.  (89)-(92)  are  depicted  in  Figs.  7a-e  for  the  cases  considered  previously  in 
Figs.  2,  4a-e,  5a-e  and  6a-e.  It  may  be  seen  from  Figs.  7a-e  that  the  functions 

and  K^*  vanish  more  rapidly  than  the  functions  and  ioK^  for  increasing  values 


of  t  *  tan8.  In  this  respect,  the  functions  K^*  and  are  comparable  to  the  func¬ 
tions  K  and  K  "  depicted  in  Figs.  6a-e.  However,  the  functions  K  *  and  K  * 
w  ri  w  rl 

depicted  in  Figs.  7a-e  and  the  functions  and  K  depicted  in  Fig9.  6a-e  are 

significantly  different  for  small  and  moderate  values  of  tan9.  More  precisely,  the 
functions  K  "  and  K  "  are  appreciably  larger  than  the  functions  K  and  ioK  ,  as 

W  rl  W 

was  already  noted,  whereas  the  functions  K  *  and  K  *  are  comparable  to  the  functions 

w  ri 

K„  and  ioK„  .  In  fact,  the  functions  K„  and  K„  are  somewhat  smaller  than  the 
W  H  W  H 

functions  K,,  and  ioK„  for  small  and  moderate  values  of  tan8. 

W  H 

Figures  5a-e  and  7a-e  show  that  the  cancellations  occurring  between  the 

waterline  and  hull  integrals  KIT  and  ioK  are  reduced  significantly  in  the  modified 

W  ri 

waterline  and  hull  integrals  K  *  and  K  *.  The  expression  for  the  Neuraann-Kelvin 

W  rl 

correction  term  K  +  +  K  “  defined  by  Eqs.  (89 )—  (91)  therefore  is  preferable  to  the 
$  4> 

expression  given  by  Eqs.  (40)  and  (41)  for  numerical  calculations. 

The  velocity  components  3<j>/3x,  3 4> /3 y  and  3$/3z  in  Eq.  (90b)  defining  the  ampli¬ 


tude  functions  A+  in  the  hull  integral  in  Eq.  (89)  can  be  expressed  in  terms  of  the 


components  <j>c  and  <j>s  of  the  velocity  vector  V<j>  along  two  unit  vectors 


(t  ,  t  ,  t  )  and  s  =  (s  ,  s  ,  s  ) 
x  y  z  x  y  z 


(93a , b ) 


tangent  to  the  hull  surface.  More  precisely,  we  have  V(j>  =  (3<(>/3n)rH-<j>tt+<J>ss  ,  which 


yields 


3<|>/3x  ■  n  3<^/3n  +  t  <)»  +  s  <j>  , 

X  X  C  X  s 

3 4> / 3 y  *  n  3<b/3n  +  t  4  +  s  (fc  , 

y  T  y  t  y  s 

3$/3z  *  n  3i/3n  +  t  6  +  s  d>  , 

z  z  t  z  s 


(94a) 

(94b) 

(94c) 


where  3<|»/3n  is  the  velocity  component  along  the  unit  outward  normal  vector  n  to  the 


hull  surface  defined  as 
n  -  (t  x  s)/  |  t  x  s  |  . 


The  unit  vectors  t  and  s  to  the  ship  hull  are  tangent  to  curves  which  approximately 
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correspond  to  waterlines  and  framelines,  respectively,  and  they  point  towards  the 
bow  and  the  keel,  respectively.  The  vectors  t  and  s  thus  are  roughly  (but  not 
necessarily  exactly)  orthogonal.  At  the  mean  free  surface,  the  vector  t  is  tangent 
to  the  top  waterline  (and  we  thus  have  t  =  0)  in  agreement  with  our  previous 
definition.  Equations  (94a-c)  yield 


n  3<b/3z-n  3 <b /3 y  =  (n  t  -n  t  )<J>  + 

y  2  y  z  z  y  t 

(n  s  -n  s  )<b  , 

y  z  z  y  s 

(96a) 

n  3<b/3x-n  9 cb /9 z  =  (n  t  -n  t  )<b  + 

2  X  Z  X  X  Z  t 

(n  s  -n  s  H  , 
z  x  x  z  s 

(96b) 

n  9<b/9y-n  9di/9x  =  (n  t  -n  t  )d>  + 

x  y  x  y  y  x  t 

(n  s  -n  s  )d> 
x  y  y  x  s 

(96c) 

By  using  Eqs.  (96a-c)  we  may  then 

express  the  function  A  defined  as 

A  =  a(n  9<|>/3z-n  3 4> /3 y )  +  8  (n  9<j)/9x-n  9<()/9z)  +  y(n  3 4> / 9 y— n  dip/dx) 
y  z  z  x  x  y 

(97) 

in  the  form 

A  =  T<t>t  +  S*s  , 

(98) 

where  T  and  S  are  given  by 

T  =  (ae  +Be  +ye  )  •  nx  t  , 
x  y  z 

(99a) 

.  .  -f  ->■ 

S  =  (ae  +ge  +yez)  •  nxs  . 
x  y 

(99b) 

By  substituting  Eq.  (95)  into  Eqs.  (99a, b)  we  may  then  express  Eq.  (98)  in  the  form 


A  =  (T'4>t  -  S"<|>  )/  |  t  x  s  |  ,  (100) 

where  T'  and  S'  are  given  by 

T'  =  as  +8s  +y  s  -  (at  +8 1  +y  t  )t»s  ,  (101a) 

x  y  z  x  y  z 

S'  =  at  +Bt  +y t  -  (as  +gs  +y s  )t*s  .  (101b) 

x  y  z  x  y  z 

Equations  (97),  (100),  and  (101a, b)  finally  yield 

a(n  9<(>/8z-n  3 4> /9 y )  +  8(n  9<j>/9x-n  9<j>/9z)  +  y(n  9 4> / 9 y— n  9(f>/9x)  = 
y  z  z  x  x  y 

[(as  +Bs  +ys  )9$/9t  -  (at  +f5t  +y  t  )9<^/9sj/  |  t  x  s  |  ,  (102) 

x  y  z  x  y  z 

where  the  relations 

9<J>/9t  =  <J>t  +  t*s  4> g  and  9<f>/9s  =  4>g  +  t«s  <f>t  (103a, b) 

were  used.  Equations  (90a, b),  (102)  and  (103a, b)  show  that  the  Neumann-Kelvm 
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correction  K  ”  defined  by  Eq.  (89)  may  be  expressed  in  terras  of  the  components 
<P 

4  and  $  of  V$  along  two  unit  vectors  t  and  s  tangent  to  the  hull, 
c  s 
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Fig.  6c  The  functions  and  for  F  *  0.2  . 
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Fig.  7b  The  functions  Kw  ,  ioK^  ,  Ky*  and  KR*  for  F  -  0.15 
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CONCLUSION 


In  summary,  the  wave  potential  <|>  (£ )  at  a  point  E,  =  (£ ,  r\ ,  q  <  0)  behind  the 
stern  of  a  ship  is  defined  by  Eq.  (20),  that  is  we  have 

4>wCO  =  (2/n)  /0  exp(v2?p2)  cos(v2npt)  Im  exp(iv2£p)  [K+(t)  +  K_(t)]  dt  ,  (104) 

where  v  and  p  are  defined  by  Eqs.  (1)  and  (14b),  respectively,  and  the  wave-spectrum 


functions  K+(t)  are  expressed  in  the  form  of  Eq.  (21),  as  follows: 


K±(t)  =  KQ ’ (t )  +  K  h(t)  . 


(105) 


In  this  expression,  represent  the  zeroth-order  slender-ship  approximation  and 

+ 

K  the  Neumann-Kelvin  correction  to  the  slender-ship  approximation. 

r 

The  slender-ship  approximat ion  K  is  given  by  Eq.  (22)  or  by  the  recommended 

modified  Eq.  (34).  We  then  have 

O(t)  =  f  E  (n  2-u2)t  dl  +  u2  /  .  exp(P2z )E  t  dJ l 
0  J  w  t  x  y  Jw  ±  y 

0  0  O') 

-  iv  u  /  exp(P  z)E,  n  da  +  v  [  exp(P  z)E.  n  da  ,  (1C 

J  s  ±  z  1  b  ±  x 


(106) 


where  w  represents  a  waterline-like  curve  separating  the  hull  side  s  and  the  hull 
bottom  b.  In  Eq.  (105),  E+  are  the  trigonometric  functions  defined  by  Eq.  (19), 
that  is 


Ef(x,y;t)  =  exp[-iP2(ux  ±  vy)]  , 


(107) 


where  P,  u  and  v  are  given  by  Eqs.  (I4a,b)  and  (I5a,b). 

+ 

The  Neumann-Kelvin  correction  terms  K  are  defined  by  Eqs.  (23)  and  (24).  This 

4> 

well-known  expression  was  modified  into  the  form  given  by  Eqs.  (40)  and  (41)  via  a 
first  mathematical  transformation.  \  second  mathematical  transformation  led  to  the 
alternative  expression  given  by  Eqs.  (74)  and  (75),  which  involve  the  arbitrary 
function  C(t).  Mathematical  and  numerical  considerations  led  to  the  selection 

+ 

Cft)  =  -uv  and  to  the  recommended  expression  for  the  Neumann-Kelvin  correction  K 

9 

given  by  Eqs.  (89),  (90a, b)  and  (91).  We  then  have 


K  '  =  fw  E±  at  dSL  t  iv  jh  exp (P  z)E±  \±  da  , 


(108) 
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(109a) 


a,  =  (tdi+S(b)t  ±  uv33d>/3t  , 

—  x  t  x  s  y 

|  t  x  s  |  A+  ■  [v3tx+u(l+v2  )t^-iuvtz]  3<f>/3s 

-  f  v3  s  +  u(l+v2)s  -iuvs  ]  3<f>/3t  ,  (109b) 

L  x  y  z 

where  Eqs.  (17)  and  (102)  were  used.  In  Eq.  (109b),  t  and  s  are  unit  vectors 
tangent  to  the  hull  surface  along  curves  which  approximately  correspond  to  water lines 


and  framelines,  respectively.  The  vectors  t  and  s  point  towards  the  bow  and  the 

keel  line,  respectively.  They  are  roughly  (but  not  necessarily  exactly)  orthogonal. 

At  the  mean  free  surface,  the  vector  t  is  tangent  to  the  top  waterline  (and  we  thus 

-*•  ,  *► 

have  t  =*0).  The  components  <j>g  and  <|>t  of  V<j>  along  the  unit  tangent  vectors  s  and  t 
and  the  velocities  3<j>/3s  =  V<j>*s  and  34> /3 1  =  V<j>»t  are  related  as  follows 

3 4> / 3 s  =  $  +e<f>t  and  3$/3t  =  ,  (110a, b) 

<j,  =  (3<t>/3 s-e3<J> /3 1 )/ (1-e2  )  and  ij>t  =  (3<j>/3t-e3<t>/3s)/(l-e2  )  (111a, b) 

where  e  is  defined  as 


e  *  t*  s  • 


(112) 


The  free-surface  integral  in  Eq.  (23)  associated  with  the  nonlinear  terms  in  the 
f ree-surface  boundary  condition  has  been  ignored  in  Eq.  (108),  which  thus  corresponds 
to  the  usual  linearized  Neumann-Kelvin  approximation.  The  generalized  Neumann-Kelvin 
expression  incorporating  the  free-surface  nonlinear  term  tt(<J>)  defined  by  Eq.  (5)  is 


then  given  by 

K^±(t)  +  If  E+  ir(<^>)  dxdy 


(113) 


The  usual  expression  for  the  functions  defined  by  Eqs.  (23)  and  (24)  involves 
both  the  velocity  potential  $  and  the  velocity  vector  V<f>.  The  alternative 
matheraat ically-equivalent  expressions  for  the  functions  K^~  given  by  Eqs.  (40)  and 
(41),  Eqs.  (74)  and  (75),  and  Eqs.  (89)  and  (90a, b)  only  involve  the  velocity  vector 


V<f>,  not  the  potential  More  precisely,  these  alternative  modified  expressions  are 
defined  in  terms  of  the  velocity  components  and  <J>g  along  the  vectors  t  and  s 
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tangent  to  the  hull,  as  may  be  seen  from  Eq.  (102)  and  is  indicated  explicitly  in 
Eqs.  (109a, b).  The  alternative  modified  expressions  obtained  in  this  study 
therefore  define  the  wave  potential  behind  the  stern  of  a  ship  in  terras  of  the  speed 
and  the  size  of  the  ship,  the  form  of  its  hull  and  the  tangential  velocity  at  the 
mean  hull  surface.  These  expressions  are  directly  suitable  for  use  in  conjunction 
with  a  boundary-integral  equation  method  based  on  a  source  distribution,  or  any 
other  numerical  method  in  which  the  velocity  vector  (but  not  the  potential)  is 
determined  on  the  mean  hull  surface. 


However,  the  main  recommendation  of  the  alternative  modified  expressions  for  the 
+ 

functions  K  obtained  in  this  study  resides  in  the  fact  that  the  cancellations 
<P 

occurring  between  the  waterline  and  hull  integrals  in  the  usual  expression  (23)  are 
considerably  reduced  in  the  modified  expressions,  especially  the  recommended  expres¬ 


sion  defined  by  Eqs.  (108)  and  (109a, b).  The  sum  of  the  port  and  starboard  contribu¬ 


tions  to  the  function  K  *  K  +  K  may  be  expressed  in  the  alternative  forms 

<f  V  <P 

K  *  K  +  ioK"  +  o2K  '  (114a) 

4>  w  W  H 

Ka  -  K  *  +  K  *  ,  (114b) 

<t>  W  H 

corresponding  to  the  usual  expression  (23)  and  the  recommended  modified  expression 


(108),  respectively.  The  waterline  integrals  K  and  ioK  '  and  the  hull  integral 

w  w 

o2K  "  in  Eq.  (114a)  are  defined  by  Eqs.  (28a-c)  and  (29),  and  the  terms  K.  *  and  K  * 

H  W  rl 

correspond  to  the  waterline  and  hull  integrals  in  the  modified  expression  (108). 

The  real  and  imaginary  parts  of  the  functions  ,  ioK.  ^  ,  V*  .  kh*  and 

K  are  depicted  in  Figs.  8a-e  for  the  simple  cases  considered  previously  in  Figs.  2, 

♦ 

4a-e,  5a-e,  6a-e  and  7a-e.  The  function  K  is  appreciably  smaller  and  vanishes  much 

<P 

more  rapidly  with  increasing  values  of  tan0  than  its  components  a^' ,  ioK  '  and  K  . 

H  W  W 

Large  cancellations  therefore  occur  among  these  components  and  the  usual  expression 


(.23)  is  quite  ill  suited  for  accurate  numerical  calculations,  notably  for  evaluating 
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the  short  divergent  waves  in  the  wave  spectrum  corresponding  to  large  values  of 

tan0.  It  may  be  seen  from  Figs.  8a-e  that  the  modified  waterline  and  hull  integrals 

K  and  K.  in  expression  (108)  are  appreciably  smaller  and  vanish  much  faster  than 
W  rt 

^  2 

the  functions  ,  ioK  "  and  a  K,/  and  are  comparable  to  the  function  K  .  The 
W  W  H  <j> 

modified  expression  (108)  thus  is  considerably  better  suited  than  the  well-known 
usual  expression  (23)  for  accurate  numerical  calculations  of  the  steady  wave  spectrum 
of  a  ship. 

For  large  values  of  t  =  tan0  ,  the  major  contribution  to  the  hull  integral  K  * 

H 

stems  from  the  upper  part  of  the  hull  surface  in  the  vicinity  of  the  mean  waterline 

due  to  the  exponential  function  exp(P  z).  The  hull  integral  K  ,  and  consequently 

the  function  K  ,  may  then  be  approximated  by  a  waterline  integral  for  large  values 
<P 

of  tan0 ,  as  has  indeed  been  shown  previously  in  this  study  for  the  special  case  of  a 
wall-sided  hull.  This  asymptotic  approximation  can  be  extended  to  arbitrary  ship 
forms,  i.e.  ships  having  flare,  and  refined  by  retaining  the  first  few  terms  in  the 
asymptotic  approximation.  A  detailed  short-wave  asymptotic  analysis  will  be 
reported  elsewhere  as  it  is  important  for  evaluating  the  short  divergent  waves  of 
interest  for  applications  to  remote-sensing  of  ship  wakes. 
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